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Abstract 

We study S- and P-wave -kit and Kir distributions for the three-body production pro- 
cesses J/ip — ► u>ir + ir~, D + — > K~tt + it + , D + — > K + tt + tt~ , and — > K + tt + tt~, by 
applying the Resonance-Spectrum Expansion to the two-body subamplitudes. The results 
are compared to the corresponding signals reported by the BES, E791 and FOCUS col- 
laborations. We find that the studied S- and P-wave channels can largely account for the 
observed structures in the data up to 1 GeV. 



1 Introduction 



Experimental data on hadronic resonances have largely shifted from elastic scattering to produc- 
tion processes over the past decade. Nevertheless, the thus observed hadrons are supposed to 
be identical to the corresponding ones first seen in elastic scattering. This generally accepted 
uniqueness of resonance mass and width, irrespective of the underlying process, reflects itself in 
the universality of the pole structure of the resonance amplitude or subamplitude. However, this 
does of course not mean that a resonance is produced with a fixed probability in all processes 
where it shows up. This will depend on microscopic details of a specific production process, the 
corresponding vertex functions, and possibly also on effects from final-state interactions with 
other particles produced alongside. 

Clearly, theory has not managed to keep up with this fast experimental development. This is 
no surprise in view of the extremely difficult dynamics at the quark level involved in a three-meson 
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decay, not to speak of decays to even more mesons. Nevertheless, it would already be significant 
progress if one could use an amplitude ansatz with the correct resonance-pole structure, as derived 
from a demonstrably working model for elastic meson-meson scattering and meson spectroscopy. 

In this spirit, we propose to start from the so-called Resonance-Spectrum Expansion (RSE) 
in order to obtain such an ansatz. For that purpose, we construct here the amplitude for the 
production of a pair of interacting mesons out of a quark- ant iquark pair, using a previously 
developed coupled-channel model [1], and its reformulation [2] in terms of an expansion in the 
radial confinement spectrum of the qq pair, the RSE. We assume thereby that the decay processes 
J ftp — > LJ7r + 7r _ , D + — > K + -k + ti~ and — > K + tt + tt~ take place in cascades of two subsequent qq 
pair-creation processes, triggered by a weak process and resulting in one spectator meson plus a 
pair of rescattered mesons. For the decay D + — > K~tt + it + , we assume that one of the 7r + mesons 
is formed in the weak process, whereas the remaining K~ir + pair stems from qq pair creation. In 
either case, the spectator meson is supposed to be formed in the first step of the cascade process, 
alongside a qq pair. The latter pair is supposed to couple to the pair of rescattered mesons. 

Accordingly, we substitute the three-meson vertex for the final process by the transition 
potential of Ref. [1]. This interaction couples qq to a pair of mesons. This pair is supposed to 
interact via quark loops, which is described by the RSE ladder sum. The quark loops, representing 
the dynamics of the two interacting mesons, can be written in a closed analytic form [1,2], and, 
consequently, the full production amplitude as well. This is not the case for the initial three- 
meson vertex — either weak or strong — which is here just represented by a complex coupling 
constant. 

The RSE has the advantage that all possible resonances for a given set of quantum numbers are 
contained in one expression. Therefore, the overlapping of resonances and their possible mixing [3] 
is automatically taken care of. Also, by employing spherical Bessel and Hankel functions for 
meson-meson distributions, the expressions become automatically real below threshold, thus 
curing the problem of subthreshold complex phases occurring in other methods [4]. But it has 
yet another advantage. The elastic two-meson (MM) scattering amplitude, which can be written 
in the form [2, 5] 

T = V + VG f V + • • ■ , 

where V represents the Born term of the transition potential V(qq — > MM) and Gf the two- 
meson propagator, leads for production to an expression of the form 

l + G f V + G f VGfV+--- = V _1 T. 

As a result, we obtain for production, in an analytically exact manner, just the common de- 
nominator of the elastic scattering amplitude for all coupled scattering channels. Hence, the 
pole structure of elastic scattering is, with our method, manifestly preserved in production. Of 
course, resonance line shapes can be different for elastic scattering and production, since those 
also depend on the numerators of the corresponding RSE expressions. 

The three-meson vertex for the final process, here represented by the transition potential 
of Ref. [1], results, through Bauer's formula [6], in a form factor determined by a spherical 
Bessel function, which depends on the angular momentum of a specific two-meson channel. This 
automatically yields the form factors and centrifugal barriers that other approaches have to 
provide empirically (see e.g. Refs. [7-9]). 

The initial vertices for strong cascades might be treated as well within the recoupling scheme 
of Ref. [10], but we shall leave this for later study. Here, the initial three-meson vertices are 
just represented by complex coupling constants, as in the isobar model (see e.g. Ref. [8] for a 
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practical application), to be adjusted to the data. Nevertheless, it should be mentioned that from 
the recoupling scheme of Ref. [10] one obtains values for the initial-process coupling constants 
which are very similar in magnitude to the ones used here. This indicates that the cascade 
assumption is very reasonable, but most probably that also some dynamics has to be considered 
for the initial processes, before all can be predicted by theory. In the prsent paper, it is our aim 
to demonstrate that the final production vertex can be handled in a very elegant way through 
the RSE. 

We shall compare our results to data taken by the BES [11], E791 [12] and FOCUS [13] 
collaborations. It turns out that the data can be very well described by our amplitudes. However, 
contrary to the FOCUS analysis [13], we shall find here that the background in their data is 
negligible. In Refs. [14,15] there are similar conclusions with respect to the data of E791 [12], 
LASS [16] and BESII [17]. It is very fortunate that such data exist nowadays. Twenty years ago 
it would have been impossible to verify whether results as presented here make any sense at all. 
Today, we may even discover that data are very well described. 

The organisation of this paper is as follows. In section 2, we discuss the details of the RSE 
for elastic scattering. The production amplitude for the RSE is derived in section 3, and the 
5*-wave predictions are compared to the data. In order to improve the description of the data 
via RSE amplitudes, we show that at least also P-waves must be included, which is discussed 
in section 4. In section 5 we give some more discussion of the details of the RSE amplitudes for 
elastic scattering. Section 6 contains our conclusions. 

2 The Resonance-Spectrum Expansion 

The RSE has been developed [2] for the description of meson-meson scattering resonances and 
bound states, in a non-perturbative approach that aims at unquenching the qq confinement 
spectrum. It consists of a simple analytic expression which can be straightforwardly applied to 
all possible non-exotic systems of two mesons. The RSE goes beyond simple spectroscopy, since 
it describes the scattering amplitude, not only at a resonance, but also for energies where no 
resonance exists. In contrast to models which have to rely upon numerical methods of solution, 
the RSE has the additional advantage that the pole structure of its scattering amplitude can 
be studied in great detail, owing to its closed analytic form. The expression for the amplitude 
can even be analytically continued, in an exact manner, to below the lowest threshold, where 
bound states show up as poles on the real energy axis. The RSE easily handles many coupled 
meson-meson channels, or coupled systems with different internal flavours. As such, it is an ideal 
expression for the study of scattering theory in general, i.e., the study of resonance structures 
and their relation to some of the many S-matrix singularities, as well as the concept of Riemann 
sheets and analytic continuation anywhere in the complex energy plane. 

The basic ingredients of the RSE are confinement and quark-pair creation. In its lowest-order 
approximation, a permanently confined quark-antiquark system is assumed, having a spectrum 
with an infinite number of bound states, related to the details of the confining force. We shall 
denote the energy levels of this confinement spectrum by E n i (n, £ q g = 0, 1, 2, . . .). Here , we 
assume that the confinement spectrum is given by 

Eni qq = E Wqq + 2nu , (1) 

which choice is anyhow rather immaterial for the purpose of our present study. The parameter 
E e represents, to lowest order, the ground-state mass of the £ q q-wave of the quark-antiquark 
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system, which is related to the effective flavour masses of the system. The strength of the 
confinement force is parametrised by u and gauged by the level splittings of the system. In 
experiment one cannot directly measure the quantities E^ q - and u, because of the large higher- 
order contributions. Consequently, the lowest-order system is purely hypothetical. Nevertheless, 
one can obtain some order-of-magnitude insight by examining mesonic spectra with more than 
one experimentally known recurrency. 

From the J PC = 1 charmonium and bottomonium states, one may conclude that the 
average level splitting is of the order of 380 MeV, leading to uj = 0.19 MeV, independent of 
flavour. The latter property is compatible with the flavour blindness of QCD, confirmed by 
experiment [18]. Indeed, the level splittings of the positive-parity f 2 mesons seem to confirm 
that flavour independence can be extended to light quarks [19]. From the ground states of the 
recurrencies one may then extract the order of magnitude of the effective quark masses, e.g. 
2m c = m(J/ip) = 3.1 GeV (in the RSE [20] we find for twice the effective charm mass the 
value 3.124 GeV), or 2m u/d = m(p) = 0.77 GeV (0.812 GeV in the RSE). For the choice (1) of 
confinement force, we determine E oiq - = m q + m q -+(1.5+£ qq -)u. Once the effective flavour masses 
and u are fixed [20], we may describe other systems, like scalar mesons and mixed flavours [21-25]. 

Through quark-pair creation the qq system is coupled to those two-meson systems which are 
allowed by quantum numbers. In principle, many different two-meson channels can couple to 
one specific quark- ant iquark system. Here, since we study the properties of the channel lowest 
in mass, we will strip the RSE of all other possible two-meson channels, thereby assuming that 
their influence far below their respective thresholds will be negligible. Via consecutive quark-pair 
creation and annihilation, a qq pair may also couple to pairs of different flavour, for instance 
uu(dd) <-> ss. Here, as we study pion-pion scattering near threshold, we shall assume that the 
coupling of a light pair of flavours to strange-antistrange can be ignored. 

The intensity of quark-pair creation is in the RSE parametrised by the flavour-independent 
parameter A. In principle, it has to be adjusted to the data. However, one may get an idea of the 
right order of magnitude by the following reasoning. For small values of A, one may determine 
the width of the ground-state resonance in the one-channel case (pion-pion here) by [26] 

sin 2 (y El -AmA 

r « A 2 £ . - , (2) 

y El -Ami 

where a represents the average distance at which light quark pairs are created, and which can 
also be defined in a flavour- independent fashion [1]. Note that this picture of a relatively well- 
defined interquark distance at which pair creation is favoured qualitatively agrees with very 
recent unquenched lattice simulations of string breaking [27]. Even our values of a, at least for 
heavy quarks, are roughly in agreement with the latter lattice predictions. For light quarks, a 
is about 0.6 fm, as we will see later on. If we take, for example, the /o(1370) resonance width 
of 0.2-0.5 GeV, then we obtain A ~ 0.60-0.95. This is of the order of 1, which would not allow 
the approximation (2). Nevertheless, we actually employ here a value for A which is of the same 
order of magnitude as our estimate (see caption of Fig. 1). 
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3 The RSE production amplitude for the 1 = tttt S- 

wave 

The RSE amplitude suitable for our purposes results from the ladder sum in quark-pair creation 
[2], and has for tttt elastic scattering the form (V and Gf represent respectively the quark-loop 
box and the two-pion propagator.) 

T(p,p) = (p\{V + VG f V+...}\p' 



2X 2 iipaj e (pa)j e (p'a) 



9nt 



■y oo 
I 

1=0 



2 Yi(2i+m(p-p') 



n=0 E nl ~ E(p) 



1 — 2i\ 2 fipa ji(pa) hf\pa) ^ 



n 2 



(3) 



n=0 E nl - E(p) 



where E n t is defined in Eq. (1), and where 



E(p) = 2^/p 2 + m 2 n and \i = n{E) = -E 



(4) 



The factors g n e are remnants [10] of the quark- ant iquark distributions associated with the con- 
finement spectrum (1). The amplitude (3) satisfies the unitarity condition |1 + 2iT\ = 1 for all 
energies E > 2771^, as can be easily verified. 
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Figure 1: Ke (red curve) and Qw(io) (blue curve) for the S-w&ve term of the RSE of Eq. (3), 
compared with the results of Caprini, Colangelo and Leutwyler (CCL) [28] (•). The RSE parameters 
are A = 1.29, a = 2.90 GeV -1 , £bi = 1-30 GeV and ui = 0.19 GeV. The recoupling parameters are here 
given by g 2 n0 = (n + 1) x 4" n . 

In Fig. 1 we compare the RSE amplitude with the dispersion-relation result of Caprini, Colan- 
gelo and Leutwyler [28]. We do not distinguish here between neutral and charged pions, and take 
the pion mass equal to the one of the charged pions. At threshold we find for the RSE 



t° (E = 2m n ) = 3?e (t° ) (E = 2m v ) 
which may be compared to data at 0.220 ± 0.005 [29]. 



0.212 



(5) 
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In Fig. 1 one observes from the behaviour of both 3fte and Sm^) that the RSE clearly 
describes a resonant structure, i.e., the /o(600) alias a meson. However, for energies above 600 
MeV, the RSE prediction does not follow the data. The main reason is the absence of a coupling 
to ss in the qq sector, as well as the KK channel (see e.g. Ref. [30]). However, for energies below 
400 MeV the agreement with the data is excellent. For more details, see Ref. [31]. 

For production processes, we need to consider the amplitude [32-34] 



A (p) = (nn \-r-k) 



d s k (nn \V\ k 



IP 



k \G f V\ p) + (k \G f VG f V\ $) + ...} 



d 6 k(nn\V\k)(k \l + G f T\p 



d 3 k(nn\V\k)(k |V _1 T|p 



(6) 



where V represents the transition potential for qq — > tttt, which, in the model of Eq. (3), is 
represented by a local delta-shell potential [2], and where V = V T [H nn — E]~ l V represents the 
quark-loop box. 

Hence, for the S"-wave meson-meson distributions in decay processes M — > spectator+meson+ 
meson, assuming that indeed final-state interactions with the spectator can be neglected, we 
obtain 



A°(E) oc A 



sin (pa) 



pa 



1 - 2A 2 /x 



sin (pa) 



pa 



,ipa 



co 77, 



E 

n=0 



1)4" 



E n — E 



(7) 



In Fig. 2 we compare the 1 = S'-wave irir amplitudes resulting from expression (7), with the 
amplitudes observed by the BES collaboration [11]. The parameters for the theoretical curves 
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Figure 2: The RSE 5-wave 1 = tttt distribution (arbitrary units), corrected for phase space (red 
curves), compared to the tttt S-wave signal of Ref. [11] extracted from J/ip — > lotttt. In (a) we compare 
the RSE to the a signal extracted from the data by BES (blue histogram), while in (b) it is compared 
to the full BES tttt signal. 

of Figs. 2 are identical to the parameters employed for the theoretical result of Fig. 1. The 
experimental data shown in Fig. 2a are identified by the BES collaboration [28] as the resonant 
1 = S-wave component, i.e., the a meson. The agreement between the RSE and the BES a data 
is reasonable, but not very good, at least for the parameter set used to reproduce the elastic 1 = 
S- wave 7i + ir~ data of Ref. [28]. For low energies, say below 600 MeV, this slight discrepancy can 
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be understood from the lower-lying a pole and somewhat larger 1 = scattering length in the 
CCL parametrisation, as compared to the BES analysis [35]. Another reason, relevant at higher 
energies, is that the RSE does not isolate resonant parts from so-called non-resonant parts of 
the scattering data. Hence, instead of forcing the RSE to agree with the a data by changing the 
parameters, we better compare the RSE to the full data. In Fig. 2b a comparison is shown of 
the RSE formula (7) with the full tttt data, obtained in Ref. [11] from J/ip — > ujtttt. 

A few words of caution are in place here. As already indicated above, the BES collaboration 
found a higher-lying sigma pole in their analysis than CCL, namely at (541 ± 39) — z(252 ± 42) 
MeV [11], to be compared to the reported [28] CCL value of (441 ± 4) - z(272 ± 6) MeV. Hence, 
the theoretical curve of Fig. 2a should a priori not follow the BES data very closely. Nevertheless, 
our prediction also seems to agree quite well, up to about 600 MeV, with the full BES mi signal in 
Fig. 2b. However, BES reported contamination of the data from processes of a different origin, so 
that no definite conclusions on pole positions and line shapes can be drawn from this comparison. 
For all that, it seems safe to conclude that the RSE is in fact capable of reasonably describing 
both the elastic scattering data of Ref. [28] and the production data of BES [11], with the same 
set of parameters, contrary to what was suggested by C. Hanhart in Ref. [36]. 

We observe that the RSE prediction follows the data [11] closely up to about 0.6 GeV, which 
is the limit of validity of the approximation (3), which neglects ss and KK contributions. The 
main differences between the RSE and the BES data stem from the D-wave ^(1270) resonance, 
and also from a small signal due to the Jo (980) resonance, both not contemplated in formula (7). 
It should furthermore be noticed that the pole structure of the amplitude in Eq. (7) is exactly 
the same as of the S-w&ve term in the elastic scattering amplitude (3). For the parameter set 
given in the caption of Fig. 1, we find the a pole at 462 — i207 MeV. 

Turning now to the kaon-pion system, in Fig. 3 we compare the 1 = 1/2 S"-wave Kir compo- 
nent of the elastic-scattering formula (3) with the data of E. Estabrooks et al. [37] as well as the 
LASS data [16], and the same component of the production formula (7) with the D + — > tc + K^tt + 
data of the E791 collaboration [12]. We obtain good agreement between the RSE and data, apart 
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Figure 3: The 1 = 1/2 S-wave Kir distribution, for elastic scattering (a) and production (b, arbitrary 
units). In (a) the prediction of Eq. (3) (red curves) is compared to the available scattering data [37] (0) 
and [16] (•). In (b) the result of Eq. (7) (red curve) is compared to the E791 analysis of D + — > K~tt + tt + 
[12]. The green curve of (a) is used in section 4. 

from a signal which is mainly due to the P-wave resonance i\"*(892). We find a contribution of 
about 41% from the ^(800) (alias kappa) resonance to the total signal. The k pole comes out 
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at 772 - z281 MeV [38]. 

Now, although the RSE amplitude (6) seems to agree well with the BES (and also the E791) 
data, it should be noticed that the experimental signal stems from a three-body decay process, 
whereas the theoretical curve in Fig. 2b concerns a two-body amplitude. So it is not obvious 
whether the comparison makes sense at all. In the following, we shall study this issue in a bit 
more detail, viz. for the 3-body decay processes 

D + -> K+tt+tt' and D+ -> K + ir + ir~ . (8) 

Here, we assume that, after the initial flavour content of the charmed mesons has been 
transformed by weak processes into us, the three-body decays under consideration take place 
as a cascade of two OZI-allowed [39] two-body processes. As a consequence, one of the three 
decay products may be considered a spectator. In this approximation, we may write the total 
amplitude in the form 

Atot = (9) 

A ([I = 0, S wave] K + ^j + A ([ I — 1/2, S wave] 7r + ) + higher waves and different isospins . 

At first, we shall neglect the higher waves, as well as possible different isospin combinations, and 
take for the total amplitude of the above decay processes just the sum of the first two terms in 
Eq. (9), i.e., terms of the form (6). In order to evaluate such an expression, we generate 10 6 
events for the allowed three-body phase space. The resulting Dalitz plot [40] is projected onto 
the A + 7r~ and 7r + 7r~ 2-body channels. In Fig. 4 we compare the RSE results to the FOCUS 
data [13] for the decay processes given in Eq. (8). 
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Figure 4: Comparison of the Dalitz-plot projections for the RSE two-particle S'-wave contributions to 
D + — > K + tt + it~ and Z)+ — ► Ar + 7r + 7r _ , and the results of the FOCUS experiment [13]. 



One difference between the RSE signals shown in Figs. 2, 3 on the one hand, and in Fig. 4 
on the other hand, is striking, namely that in the latter (production) case the amplitudes do not 
vanish anywhere except for the boundaries of allowed phase space. The reason appears to be 
simple. Namely, each of the two S'-wave amplitudes of formula (9) has its own zeros associated 
with it, which lie at the energies of the corresponding confinement spectra, and are different 
for different flavour combinations (see formula (1)). Hence, the sum will in general not vanish, 
particularly not for a Dalitz-plot projection at a constant invariant mass in one of the two S'- 
wave amplitudes, integrated over all phase-space-allowed invariant masses for the other S-wave 
amplitude. 

The data have been collected in equal bins of invariant mass squared by the FOCUS collab- 
oration, and so are the RSE results. However, we prefer to present the results as a function of 
mass, rather than mass squared. Instead of manipulating the data for that purpose, we let the 
bin size vary with mass in Fig. 4. The RSE results for the processes (8) have been normalised by 
hand, in order to be comparable to the FOCUS data [13]. Our choice of normalisation suggests 
that the dominant modes are given by two-body S-waves. The remaining structures can visibly 
be associated with the /o(980) (dominantly ss) component of isoscalar ++ states (here, we only 
consider the nn component), as well as P- and D-wave contributions, like the p(770) and the 
/ 2 (1270) in vr+vr-, and the if* (892), AT* (1410), A'* (1430) and A'* (1680) in AT+vr", besides their 
interference patterns with the S-wave signals. 
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The RSE amplitude (9), which consists of the sum of two amplitudes of the form (7), has the 
same pole structure as the individual S'-wave isoscalar -kit and isodoublet Ktt elastic-scattering 
amplitudes. So one can only determine the pole positions in the production amplitude once the 
model parameters have been fixed. Moreover, one cannot just remove a pole from this amplitude, 
in order to verify whether or not it contributes significantly to the experimental signal. Whereas, 
in most cases, one can relate resonant structures in elastic scattering to each of the poles in 
the scattering amplitude, here such a relation is not evident at all, since most of the poles [38] 
in the and Kit S-wave amplitudes do not show up as peaked structures in Fig. 4 for the 
total RSE amplitude (9). In fact, most of the theoretical S-wave signal has the aspect of a 
flat non-resonating background, except for the clear a and k signals at 470 MeV in 7T + ir~ and 
800 MeV in K + n~, respectively. Nevertheless, the amplitudes (7) contain all poles of elastic 
S'-wave scattering. Furthermore, it is not easy to find out whether the theoretical structure in 
Kir around 1.6 GeV stems from the lower pole, which corresponds to the _Ko(1430) resonance in 
elastic scattering, or from a higher pole. 



4 Study of the P wave contributions 

Although the RSE results shown in Fig. 4 agree well with experiment [13], it must be acknowl- 
edged that in the Ktt Dalitz-plot projections at low energies, which is precisely the mass region 
where the RSE is supposed to work well, the predicted amplitudes for — > K + tt + tt~ and 
D + — > K + tt + tt~ are larger respectively much larger than those reported in Ref. [13]. Hence, the 
legitime question arises whether the interferences of the tttt and Kir S-waves with higher waves 
are capable of correcting these discrepancies, as claimed in the previous section. Thereto, we 
embark on an extension of expression (9) to P- waves, i.e., 



.4 



tot 



A MI = 0, S wave] K + ) + A ( [I = 1/2, S wave 



7T 



(10) 



A ([I = 0, P wave] K + ) + A ( [I = 1/2, P wave] vr 



The generalisation [2] of the RSE two-particle one-channel amplitude (7) to higher waves is given 
by 



A e (E) oc X qq i e je(pa qi! ) 



2iX 2 q -npa qq j e (pa qq ) hfp (pa qq ) 



9nl 



n=0 E n — E 



where [10] g^ = (n + 1)4 n and g\ x = |(2ra + 3)4 n , and where a qq and X qq can be related to the 
flavour-independent 3 P coupling A and average transition distance a by 



A 



qq 



uj (m q + m g 
m q m q 



A 



:i2i 



The procedure to calculate the total amplitude for each of the possible invariant-mass com- 
binations of tttt and Kir, is as follows. First, we adjust the only "free" parameters, viz. A and a, 
of each of the two-particle P waves to the elastic scattering data, such that the p and .ft"* (892) 
resonances are well reproduced. Then, we adjust the fraction with which each of the two P-wave 
amplitudes contributes to the total amplitude (10). Apart from an overall phase of 90° between 
the P- and S-waves, there is no need to introduce additional relative phases. Hence, in total, the 
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Figure 5: Comparison of the Dalitz-plot projections for the RSE two-particle S- and P-wave contri- 
butions to D + —¥ K + tt +/ k~ and Df — > K + ir + n~, and the results of the FOCUS experiment [13]. The 
total number of events under the theoretical curves is 10 percent less than in experiment for D + decay, 
and 20 percent for Df, in order to account for the signal which is missing in the theoretical curves at 
higher energies. 
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relative contributions of the four terms in expression (10) can be parametrised by 3 real numbers. 

It is clear from the resulting amplitudes shown in Fig. 5, that the interference between S- 
and P-waves can perfectly account for the discrepancies observed in Fig. 4. In particular, the 
k resonance seems to have disappeared from the D + — > K + tc + tt~ signal, in agreement with 
experiment. Nevertheless, the RSE amplitude for D + — > K + 7T + ir~ does contain the k pole, 
namely at 758 — z279 MeV. Consequently, its apparent absense from the signal does not mean 
that it does not contribute to the amplitude. 

In Table 1 we show the relative intensities of the four contributions of Eq. (10). 



relative 


(K + TT-) S TT + 


(tt + it-) p K + 


(K + n~) p n + 


ratios 


(7T + 7r-) s K + 




(it + ti-) s K + 


for D + decay 


1.5 


1.7 


1.4 


for Df decay 


2.5 


2.0 


0.7 



Table 1: The relative intensities of the four amplitudes of Eq. (10). 

In Fig. 6, we also find improvement at low energies with respect to Fig. 3b, by including the 
P-wave amplitude for the process D + — > K~7T + ir + , measured by the E791 collaboration [12]. 
The data of Fig. 6 were collected by E791 by selecting those K~n + invariant masses which are 
lowest in mass for the two possible K~n + combinations in the decay process D + — > K~7T + 7r + . 
In our Monte-Carlo event generation we perform the same selection procedure for the RSE. 




0.8 1.0 1.2 

™-K-7r+(GeV) 



Figure 6: Comparison of the Dalitz-plot projections for the RSE two-particle S- and P-wave contri- 
butions to D + — > K~ir + 7r + , and the results of the E791 experiment [12]. The relative intensity of the 
P-wave contribution is 0.35 with respect to the 5-wave contribution, and has a relative phase of 180° 
for the theoretical curve. The RSE parameters are taken from Table 2. 

We recall that the production processes D + K + tt + tt~ and D + — > K~tt + tt + are different in 
nature. This difference reflects itself, in particular, in different relative phases, viz. 90° for the 
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former process and 180° for the latter. Hence, the line shapes of Figs. 5 and 6 are very different. 
Nevertheless, as in the D + — > K + ir + 7r~ cross section, also in the D + — > K~ix + ix + signal there 
is no clear presence of the k enhancement, although the RSE amplitude contains the k pole at 
758 - i279 MeV. 

5 Further discussion of the RSE 

Formulae (7) and (11) are very incomplete expressions for the amplitudes of meson-meson pro- 
duction. First of all, they are based on a set of dynamical equations in which the interaction 
describing the transition from quark-antiquark to meson-meson, via the 3 Po mechanism, is mod- 
elled through a one-delta-shell potential with radius a. The advantage of such an approach is that 
the confinement spectrum (1) enters the expressions explicitly. However, it has the consequence 
that a must be adjusted for each separate configuration, thus losing universality. Secondly, by 
reducing to one the number of possible channels to which a qq state is allowed to couple, one 
must adjust also A, in order to get the elastic-scattering poles in the correct positions. So one 
also loses the universality of the 3 Pq coupling. On the other hand, the expressions for the am- 
plitudes become much less transparent when all ingredients of flavour independence and channel 
completeness are built in. In the present paper, our principal aim is to show that the RSE re- 
produces the basic properties of meson-meson production, and not to carry out a detailed data 
analysis. The price we have to pay is that a and A are not completely constant for the four 
amplitudes in question. In Table 2 we give the values used for the curves in Fig. 5. 



model 
parameters 




(K + tt-) s it + 


(lT + lT-) p K + 


(K + TX~) p Tl + 


\a 


1.00 
1.000 


1.09 
0.654 


1.12 
0.960 


1.01 
0.995 



Table 2: The values of a and A in each of the four amplitudes of Eq. (10). 

The flavour-independent coupling constant A represents the intensity of qq pair creation and 
annihilation, which is considered a constant in the RSE in case all meson-meson channels allowed 
by quantum numbers are taken into account. In practice, the latter is an impossible task, of 
course. However, omitting meson-meson channels has an influence on the resulting resonance pole 
positions. Thus, the coupling constants of the remaining channels act as effective parameters, 
which must be adjusted in order to reproduce the elastic scattering data. At low energies, there 
is a negligible effect on pole positions and resonance shapes from meson-meson channels with 
much higher thresholds. Accordingly, their omission has little influence on A. So, for energies 
below 1 GeV, it is a reasonable approximation to restrict the allowed meson-meson channels to 
those containing a pair of pseudoscalar mesons, and even to only irir and K7i, for the isoscalar 
and isodoublet cases, respectively. The only situation where this philosophy does not work well 
here is for the k(800) pole. The reason is that the Kv( channel takes [41] about half of the 
total coupling, but nevertheless does not contribute significantly to the k pole position, since this 
high-lying channel is strongly damped there. As a consequence, in this case the effective coupling 
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for Kit to J p = + ns comes out some 35 percent below the average flavour-independent value 
of A. 

The effective average interaction radius a mainly controls the resonance width, which is 
relatively sensitive to the omission of higher thresholds, as can be observed from Table 2. 

The effective flavour masses m u /d = 406 MeV and m s = 508 MeV, as well as the average 
mesonic level splitting uo = 0.19 GeV, have been determined in Ref. [20]. These parameters 
generate the quenched quark- ant iquark spectrum, which should be compared to corresponding 
lattice calculations. For example, in a very recent preprint [42] the Scalar Collaboration found 
1.58±0.35 GeV for the lightest isoscalar scalar meson and 1.71±0.54 GeV for its isodoublet 
partner, which in the RSE (Ref. [21], 1986) came out (Eq. (1) for Poi = m q + m q + 5u;/2) 
at 1.29 and 1.39 GeV, respectively. However, by assuming the latter values in the quenched 
approximation, one can in the RSE even determine the lowest S"-matrix poles [21], by adjusting 
the model parameters a and A to the available elastic S- [21] and P-wave [20] scattering data. 
Moreover, as mentioned above, the very string-breaking picture underlying the RSE, and even 
the empirically found values for a have received strong support from recent unquenched lattice 
calculations [27]. Therefore, the RSE is really capable of relating experiment to QCD, via the 
lattice. 

In previous work [38], and also for the theoretical curves in red shown in Fig. 3, we used a 
value of 1.31 GeV for the ground state of the J p = + ns confinement spectrum, which is about 
80 MeV below its value following from Eq. (1), in order to account for the contribution of the Krj' 
channel. Here, we have restored this, leading to small adjustments for a and A, and to the green 
curve in Fig. 3 for the elastic RSE amplitude. The corresponding k pole thus moves from [38] 
772 - «281 MeV to 758 - i279 MeV. 

The "poles" of the confinement spectrum, which enter formula (11), are then fully determined 
by the parameters of Ref. [20], i.e., E n e qq = m q + m q + (2n + £ qq + 1.5)u;. This amounts, for the 
5"-wave amplitudes of Eq. (10) (£ qq = 1), to 1.29, 1.67, 2.05, . . . GeV for the isoscalars, and 1.39, 
1.77, 2.15, ...GeV for the isodoublets. Notice that the k pole lies about 600 MeV below the 
lowest "pole" of the input confinement spectrum (see also Fig. 1 of Ref. [43]). The reason is 
that resonance poles do not occur for E = E n in formula (11), but only when the denominator 
vanishes. In Ref. [44] pole trajectories were shown for variations in the parameter A of formula 
(11). 

For the P-wave amplitudes of Eq. (10) (£ qq = 0), we have a spectrum at 1.10, 1.48, 1.86, 
. . . GeV for the isoscalars, and 1.20, 1.58, 1.96, . . . GeV for the isodoublets. Furthermore, for the 
parameter set of Table 2, we find the a pole at 462 — i207 MeV (also see section 3), the p pole 
at 756 — 266. 5 MeV, corresponding to a central resonance peak at 772 MeV and a width of 151 
MeV, and the K* pole at 892 - i24.8 MeV, resulting in a peak at 895 MeV with a width of 52 
MeV. 

6 Conclusions 

The spectrum entering the RSE amplitudes as input is that of a confined ("quenched") quark- 
antiquark pair. In the present work, this spectrum is parametrised by the non-strange and strange 
effective quark masses m n and m s , and the strong oscillator frequency u, all taken from Ref. [20]. 
The mesonic resonances and their three-meson vertices are fully generated by the RSE, and not 
input as in related yet perturbative methods. 

The predictive power of the RSE as an analytic method to unquench the quark model has 



11 



been demonstrated before, by interrelating an enormous variety of non-exotic mesonic systems, 
such as the light scalar mesons /o(600), /o(980), J<"q(800), ao(980) and the corresponding S'-wave 
mi, Kir, t]ti scattering observables [21,30], the scalars between 1.3 and 1.5 GeV [21], vector 
and pseudoscalar mesons [20], charmonium and bottomonium [1], the -D* (2317) [22], and the 
_D*j(2860) [25]. In the present work, we have shown that also for tttt and Kir production, i.e., in 
three-body decay processes, the RSE is a very powerful method to describe and understand the 
experimental data. 
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